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We introduce a new approach to analyzing the interaction between classical and quantum systems 
that is based on a limiting procedure applied to multi-particle Schrodinger equations. The limit 
equations obtained by this procedure, which we refer to as the classical-quantum limit, govern the 
interaction between classical and quantum systems, and they possess many desirable properties 
that are inherited in the limit from the multi-particle quantum system. As an application, we 
use the classical-quantum limit equations to identify the source of the non-local signalling that is 
known to occur in the classical-quantum hybrid scheme of Hall and Reginatto. We also derive the 
first order correction to the classical-quantum limit equation to obtain a fully consistent first order 
approximation to the Schrodinger equation that should be accurate for modeling the interaction 
between particles of disparate mass in the regime where the particles with the larger masses are 
effectively classical. 
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I. INTRODUCTION 

Consistent schemes for the interaction of classical and quantum systems are of great interest from both theoretical 
and practical points of view. For example, in the Copenhagen interpretation of quantum mechanics, measuring devices 
are treated as classical objects that interact with quantum systems. Consequently, from a foundational perspective, 
having a consistent and effective method for analyzing the interaction between classical and quantum systems is 
highly desirable. From the perspective of computational chemistry, the ability to treat heavy particles classically 
and light particles quantum mechanically results in a significant reduction in computational costs over fully quantum 
mechanical treatments. For systems with large particle numbers, this separation is essential to making computations 
feasible. 

There has been much work done on the subject of the interaction of classical and quantum systems, and various 
schemes have been developed to investigate these types of interactions. In the literature these are known as classical- 
quantum hybrid schemes, and a wide range of distinct schemes have been developed for various purposes including 
modeling of chemical reactions, analyzing decoherence, and studying measurement theory 

Additionally, a number of studies of the general structure and properties of classical-quantum hybrid schemes have 
been carried out [I|, H, IflI-llTl [2bl - l^ . 

Although many classical-quantum hybrid schemes have been introduced, none are regarded as being the definitive 
solution to the problem of describing the interaction between classical and quantum systems. In this article, we 
introduce a new method to address this probleiu that is based on taking a singular limit of a fully quantum system 
in complete analogy with the more well-known classical limit. More specifically, we introduce a small dimensionless 
parameter e into a multi-particle Schrddinger equation in such a way that the parameter only appears in a subset 
of the single particle Hamiltonian operators. Using a non-local polar representation for the wave function, we show 
that the multi-particle Schrodinger equation, when expressed in terms of our non-local polar decomposition, admits 
a formal e 0 limit. The resulting limit equations define our notion of a classical-quantum limit, which as far as we 
are aware is new. The derivation of these equations is carried out in Section [Hi and we also justify our interpretation 
of these limit equations as describing the interaction between a classical and quantum system in that section. 

The limit equations derived here define a completely natural classical-quantum hybrid scheme that has many 
attractive properties. We single out two properties. First, since our limit equations arise from the taking the singular 
limit of a fully quantum mechanical system, the desirable properties that hold true for the quantum system remain 
true for limit equations for any quantities with a well defined limit. In particular, this implies that the limit equations 
that arise from a non-interacting quantum system will satisfy the non-signalling property that is known to be satisfied 
by non-interacting quantum systems. This fact forms the basis of our first application where in Section Hill we use our 
classical-quantum limit equations to clarify the origin of the non-local signalling [T^ that was observed to occur in 
the classical-quantum hybrid scheme of Hall and Reginatto [l3, Second, again due to the limiting procedure, 

it can be expected that solutions of the limit equations will approximate corresponding fully quantum mechanical 
solutions, at least on suitable time scales as is known to be true in the classical limit Q. From this perspective, our 
limit equations capture the interactions between classical and quantum systems in a manner that is fully consistent 
with quantum theory. Moreover, we can view our limit equations as the zeroth order equations satisfied by an 
e expansion of e-dependent solutions to the fully quantum mechanical equations. Using a standard perturbation 
expansion method, we derive in Section fiVI the equations satisfied by the first order expansion term. In this way, it is 
possible to derive a fully consistent perturbative approximation scheme for multi-particle quantum systems where the 
base approximation consists of an interacting classical-quantum system. We expect that this scheme will be useful in 
computational chemistry for modelling the interaction of heavy particles that are treated in the first approximation 
classically with light particles that are treated quantum mechanically. 


II. A CLASSICAL-QUANTUM LIMIT 


Our derivation of the classical-quantum limit begins with a 2-particle Schrodinger equation of the form 

ihdtilj =+ UiIj - + Vip, (II.l) 

where U = U (x) and V = V(x, y) are time-independent potentials and we are using dx and dy to denote partial 
derivatives, i.e. 


dx(t> 


d(j) 

dx 


and 



for any function (j) = 4>{x, y). 
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Remark II. 1. In the following, the partial derivatives will act on functions that directly follow them. For example, 

OX 

If we want the partial derivative to act on both functions, then we will employ brackets and write 




dx 


Next, we introduce time and length scales given by T, Li and L 2 , and define the dimensionless variables 


t = 


t 


X = 


y = 


y 


T’ Li’ " L2 

= ^{Ti,LiX,L2y), 
TMo 


and 


V{x,y) = —V{Lix,L 2 y). 

Writing the Schrodinger equation inu in terms of these dimensionless variables, we see that 'll) satisfies 

hT r Ml -- r hT 


ids'll; = -- 


2MiLl 

We introduce a small dimensionless parameter e via 


Ml, 

M 2 


^2'*/' + 777^-72 


2M2LI y- 


■V'll). 


e = 


M 2 


(II.2) 


(II.3) 


Defining the dimensionless masses 


mi = 


M 2 LI 


and 1712 = 


12^2 


hT HT ' 

and dropping the tildes, we can write the dimensionless two-particle Schrodinger equation (III. 21) as 

1 


idtip" = - 


^dir + -ur - 

2mi e 




(II.4) 


where we have written ip = to emphasize that the solution depends on the small parameter e. We remark that we 
restrict ourselves to a pair of I-dimensional particles and time-independent potentials for simplicity. All of our results 
generalize to higher dimensions, more particles and time-dependent potentials in an obvious manner. 

The equation (III.4|) governs the evolution of a pair of 1-dimensional quantum particles of (dimensionless) mass wi/e 
and 7712, respectively. Assuming the length scales Li and L 2 are comparable, the ratio (r77i/e)/r772 ~ 1/e is unbounded 
as e 0 showing that the limit e 0 is an extreme mass ratio limit that represents the interaction of a heavy particle 
with a light one. Furthermore, due to the positioning of the parameter e, it is clear that the e \ 0 limit acts like a 
classical limit, but only for the first system with the configuration coordinate x. From these considerations, we expect 
that a consistent e 0 limit represents the interaction of a classical particle with a quantum one. We justify this 
expectation in the remainder of this article. 


A. A non-local polar representation 

The classical-quantum limit considered here is based on a non-local polar representation for the wave function where 
the non-locality is introduced by way of the averaging operator 


A{(p){x,y) := 


(b-a) 


(p{x, z) dz 


(II.5) 
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with a, & G M fixed constants satisfying a < b. We observe that A satisfies 

= i, (11.6) 

and hence, defines a projection operator. It is also easy to see that A satisfies the properties 

dxA = Adx and dyA = 0, (II-7) 

and that 

A(fi = (j) and A{(l)^) = (j^A^ for cj) = 4>{x) and ^ = ^{x,y). (II-8) 

Remark II.2. The properties (III. 61) - (III. 811 are the essential properties of the operator A that are needed to define our 
classical-quantum limit. Given any operator A satisfying these properties, all of the arguments used in this article 
remain valid, and hence, there exists a classical-quantum limit associated with every operator A satisfying (111.611 - 
(III. 81) . and in particular, for each choice of a, 6 G K, a < &. These different classical-quantum limits are all related 
via a bijective transformation, and consequently, different choices of A lead to different representations of the same 
classical-quantum limit', see Appendix 1X1 for the details. Thus there is no physical content in the particular choice 
of operator A satisfying (III.6|) - (III.8I1 . and throughout this article, we will employ the operator defined by (|II.5II for a 
fixed choice of constants a, 6 G M, a < 5. 

Exponentiating A yields 

exp(TA) = e^A -\- B (II-9) 

where 

.B := I - A 

is the complementary projection operator, which, by (III. 61) . (III.7I1 and (III. 81) . satisfies 

B^ = B, bA = Ab = q, A + B = 1, (II.IO) 

d,B = Bd^, dyB = dy, (II.ll) 

and 

B(j) = 0 and B{(j}^) = (pBS^ for (j) = 4>{x) and ^ = ^{x,y). (11-12) 

Letting 

v) = y) {RAx, y) := \r{x, y)\) (11.13) 

denote the standard polar representation, we introduce a non-local phase variable 9'^ via the definition 

9^ := exp(ln(e)i)0'^ + BQ\ (11.14) 

or equivalently 

0*^ = exp(- \n{e)A)9^ ^ B9\ (11.15) 

e 

Substituting (III.131) into (III. 41) . we find using (III. 151) that the pair {R'^,9A satisfies 

dtR" = — —d„;{exp{-\n{e)A)9'')dj;R'' - —^a^(exp(-ln(e)A)6('')i?'' 
mi Zmi 

-^i9„(exp(—ln(e)A)0'^)9„i?'^ — ■-^9^(exp(—ln(e)A)0'^)i?'^, (11.16) 

TO2 2m2 

dt9^ = exp(ln(e)i)| - ^ [9, (exp(-ln(e)i)0^)]" - 


(11.17) 
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which is completely equivalent to the 2-particle Schrodinger equation (III. 41) for any e > 0. We claim that this system 
has a regular e \ 0 limit. To see this, we observe, from repeated use of the relations (III.6I) - (III.12I) . that the following 
terms that appear to be singular are in fact regular: 


aj,(exp(-ln(e)i)r) = = dySO^ = dy9\ 


and 


exp(ln(e)A)-17 = - {eA + B)AU = U, 
e e 


exp(ln(e)y4)e[9x(exp(—ln(e)^)0'^)]"' = e(tA + B)( ^\^xA6^^^'^ + -dxAO^dxBO^ + [dxBO'^ 


(11.18) 

(11.19) 


= {dxAO'^y + e^A{dxBe^f + 2dxAe^dxBe + eB{dxB9^y 


( 11 . 20 ) 


Using l|11.18() . pi.191) and (111.201) along with the relations (I11.6I) - (I11.12I) . it is not difficult to verify that (I11.16|) - pi.l7l) 
can be written as 


dtW = — -(dxA9'^ + edxB9AdxR^ -—{dlA9^ + - —dy9^dyR^ - —d^9R<^, 

mi ^ ’ 2mi ^ ^ ! m2 2m2 " 


dt9^ = - 


1 


2mi - 


{dxA9'^y + €^A{dxB9^y + 2dxA9^dxB9^ + eB{dxB9^y 


- U 


+ {eA -f B) 


dlR^ 


1 


2mi R'^ 2m2 


(a,D‘) 


», 1 


2m2 R^ 


Setting 


9\ := A9^ and 9% := 


( 11 . 21 ) 


( 11 . 22 ) 


(11.23) 


we see after applying the projection operators A and B to (111.221) . and again using the relations (I11.6I) - (I11.12I) that 
(lll.21l) - (lll.22l) are equivalent to 


dtR" = -:^{9 x 9 a + edx9%)dxR'' - Tj^{dl9A + edl9g)R'' - -^dy9gdyR'' - 


dt9\ = - 


TOi 

1 

2mi . 


{dx9\f + e^A{dx9y) 


2mi 
2 


-U + eA 


m2 

e dlR"^ 1 


dt9% = B 


2m ^{px9y) ^ -f 


2mi R'^ 2m2 

e dlR^ 1 


[dyWy 


2mi R'^ 2m2 


{dy0%y 


2 , 1 dyR^ y 

2m2 R^ 

1 dlR^ y 


2m2 R'^ 


(11.24) 

(11.25) 

(11.26) 


It is worth noting that we can treat 9\ — {t, X, y) as independent variables provided that the 

initial data for ds at t = 0 is chosen to satisfy the constraint A9'^g\t=o = 0. This is possible because it follows easily 
from (111.261) that the constraint A9^g = 0 propagates, i.e. satisfies dtA9% = 0, from which it follows that the function 
9'^ defined by 9'^{t,x,y) = 9A{t,x) + 9g{t,x,y) verifies the relations 9% = A9‘^ and 9g = B9^ in accordance with 

(IIM . 

To summarize, the two systems (lll.21l) - (lll.22l) and (lll.24l) - (lll.26l) represent equivalent formulations of the 2-particle 
Schrodinger equation (111.41) under the transformation 

^ (11.27) 

where 


9^ = 9\+9%, (11.28) 

and 9\ and 9% are given in terms of 9'^ by (111.231) . Moreover, it is clear from (111.271) . that if yn is the characteristic 
function of a subset U C i.e. Xn(a:,J/) = 1 if {x,y) G U and 0 otherwise, then the action of the operator XQ{x,y) 
on the wave function is equivalent to the following action on the pair (0*^, R'^): 

XQ.{x,y){9\R^){x,y) = {9''{x,y),xn{x,y)RAx,y)). 

From this, we conclude that position measurements are represented by e-independent maps in our non-local polar 
representation. 
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B. The limit equations 


Taking the formal e \ 0 limit (III.21I) - (IIL22I) and (IIL24I) - (III.26I) . we obtain the two equivalent forms of the limit 
equations given by 


dtR = -—d^{A0)d^R - ^dliAe)R - —dy0dyR - ^dl0R, 

2mi m 2 2m2 


dt0 = - 


mi 
1 

2mi - 


{d^A0f + 2d^A0d^B0 


-U + B 


1 


2m: 


\dy0) 


2 1 dlR ■ 

I -V 


2m2 R 


(11.29) 

(11.30) 


and 


dtR = -—d^Ad^R - -^dlOAR - —dyOedyR - ^dl0BR, 

mi 2mi m 2 2 m 2 " 

dt^A = ~ 


dt0B = -—d^0Ad,0B + B 
mi 


1 


2m2 


(dyOB) 


2 , 1 d^R ■ 

I-«- y 


2m2 R 


(11.31) 

(11.32) 

(11.33) 


respectively. In the second formulation given by (III.31ll - (III.33|) . we can, as above, treat 0a = 0A{t,x) and 0b = 
0B(t,x, y) as independent variables provided that we choose initial data for Ob satisfying A0B\t=o = 0. 

These limit equations define our notion of a classical-quantum limit. The classical nature of the limit system is 
clear from equation (III.321) . which is of the form of a Hamilton-Jacobi equation from classical mechanics. We note 
that it only depends on the configuration coordinate x of the first system, which allows us to interpret the first system 
as being classical and x as a classical configuration coordinate. Further, we observe that the first equation (III.31|) is 
the conservation law for the probability density 


p = R^ 

that is common to both quantum and classical systems. This can be seen by introducing the velocity field 

V = (z;i,t>2) := (—d^0A, —dy0B\ 

\mi m 2 ) 

which allows us to write (|II.31I) as the conservation law 

dtp A div(vp) = 0. (11.34) 

It is also apparent that equation (III.331) is closely related to the quantum mechanical equation satisfied by the phase 
function with the main differences being the appearance of the complementary projection operator B and the velocity 
coupling term —i^dxdAdx^B = —Vxdx0B- Without these terms, (III.331) would be identical to the quantum mechanical 
equation for the phase function, i.e. a Hamilton-Jacobi equation with a quantum potential. As we show below in 
Section imi by using Lagrangian coordinates adapted to the velocity field Vx along with a change of variables, 
equations (III.311) and (III.331) can be cast in the form of an exact Schrodinger equation, albeit with a modified time- 
dependant potential. This justifies the interpretations of the second system as being quantum mechanical. 

An important point regarding these limit equations is that we can expect that solutions to these equations to 
approximate solutions to order e of the fully quantum system on time intervals [0,T] for which the solutions to the 
limit equations remain smooth. This is because the limit equations are derived from taking the e ^ 0 of the e- 
dependent 2-particle Schrodinger equation (Hn, and for related limits, such as the classical limit, there exists many 
rigorous mathematical results of this character available in the literature; for example, see Q and reference therein. 
Thus, we expect that if (i?'^,6i^,0^) is a solution of the full quantum mechanical system, then 

R^ = R 0(e), 0^ = 0A + 0(e), and = 0b + 0(e) 

on some hxed time interval [0,T] where {R,0a,0b) is a solution of the limit equations (III.31I) - (III.33I) . In this way, a 
solution (R, 0a,0b) of the limit equations can be interpreted as the zeroth order expansion term in an e expansion of 
a quantum mechanical solution of the e-dependent 2-particle Schrodinger equation dull). From this starting point, 
higher order expansion terms in e can be derived in a standard fashion. We derive the equations satished by the first 
order correction in Section llVI 
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From similar results for the classical limit (see i), we can expect that time interval [0,T] on which solutions to 
the limit equations will approximate the fully quantum system will be determined by the time interval for which the 
limit equations admit smooth solutions. As we establish below in Section fllDl this time interval is determined by the 
time of existence for the classical Hamiltonian flow associated to the Hamilton-Jacobi equation (III.321) and the time of 
existence of smooth solutions to (III.321) . which is well known to be governed by the formation of caustics. Presumably, 
in line with known results for the classical limit, see again Q, it should be possible using a more complicated multi¬ 
phase construction to extend the time of the validity of the approximation to arbitrary time intervals. We will not 
consider this aspect further here and instead leave it for future consideration. 


C. No quantum backreaction 

In the setting of quantum-classical hybrid schemes, the notion of quantum backreaction refers to the effect of the 
quantum subsystem on the classical one. Here, we show that the classical-quantum limit equations derived above do 
not allow for quantum backreaction. This is in concordance with the theoretical studies carried out in [^, [26l - [^ . It 
also makes complete sense from the origin of the limit equations as an extreme mass ratio limit. In such a limit, one 
does not expect the light particles to affect the heavy ones at the “zeroth order”. However, it is worth noting that 
for general choices of the potential V = V{x,y), the classical equations do have an effect on the quantum ones, and 
so, in general, there is a non-trivial interaction between the two systems. 

It is important to note that it is still possible to augment our limit equations with interactions that do allow for a 
quantum backreaction on the classical subsystem. Indeed, as it will be shown in Section IHII by adding appropriate 
interaction terms to the limit equations, we can obtain the Hall-Reginatto hybrid scheme equations which do allow 
for a quantum backreaction [3 [mu. Of course, this means leaving the realm of quantum mechanics and limits 
thereof into the regime of classical-quantum hybrids, which are theories distinct from quantum mechanics as we know 
it. On the other hand, a quantum backreaction of the quantum system on the classical one can be obtained within 
the framework of quantum mechanics by considering the first order correction to the limit equations that are derived 
in Section lYl From our point of view, this is a more fundamental and promising direction for understanding the 
“true” interaction between nearly classical systems and fully quantum ones. 

In order to establish that our classical-quantum limit equations do not generate a quantum backreaction, we consider 
the marginal density 

= / p{t,x,y)dy, 

Jr 

which contains all of the statistical information relating to the classical system. It follows from integrating (III.341) 
together with the fact that vi is independent of y that pi satisfies 

dtpi + dxivipi) = 0 . 

Noting that this equation and pi.321) form a closed system of evolution equations for (OatPi), it follows immediately 
that the quantum subsystem does not have any effect on the classical system and establishes the absence of a quantum 
backreaction. 


D. A local formulation 


Due to the appearance of the averaging operator in the limit equations, they are manifestly nonlocal, and conse¬ 
quently, difficult to interpret, and from a computational point of view, are more expensive to solve compared to local 
systems. Fortunately, solutions to the limit equations can be obtained from solutions of the following local system: 


dtR = 
dteA = 


dtOs = 


-—djAd.R - ^BIOaR - —dyOedyR - -^dlHR, 

mi 2mi TO2 2m2 " 


— ^O^OaOxOb - 

mi 


1 

2m2 


{dyOBy + 


1 d^R 
2m2 R 


-V. 


(11.35) 

(11.36) 

(11.37) 


Indeed, given a solution 


(0A = 0A(t,x),dB = 0B(t,x,y),R = R(t,x,y)) 

















of pi.35|) - (IIL37|) . a straightforward computation using the relations (III. 101) - pi. 121) establishes that 

{eA,eB ■■=Mb,r) 


(11.38) 


is a solution of (lll.31l) - (lll.33|) . 

We can further simplify the local system (lll.35l) - (lll.37l) by introducing a Lagrangian coordinate for the classical 
subsystem that is adapted to the vector field 


vi{t,x) = — dx0A{t,x), 

mi 


(11.39) 


which determines the motion of the classical particle. The Lagrangian transformation is defined using the flow F(t, x) 
of vi(t,x), which is uniquely determined by the initial value problem 


dtF{t,x) = vi{t,F{t,x)), 
F(0, x) = X. 

We note that the initial condition (111.411) implies that 

dxF{0,x) = 1, 

and this, in turn, guarantees that the condition 

dxF{t, x) > 0 


(11.40) 

(11.41) 


(11.42) 


holds on some time interval [0,T]. The condition pi.42l) is precisely the necessary condition for the Lagrangian 
coordinate transformation to be invertible; for the remainder of this section, we will assume that (I1L42|) is satisfied. 
Defining 


aA(t,x) := 9A{t,F{t,x)), 
aB{t,x,y) := 9B{t, F{t, x),y) 


and 


n{t,x,y) := R{t, F{t, x),y)y^d^F{t, x), 
it follows directly from (I1L39L pL40l) and the chain rule that 

dtaAit,x) = dt9A{t,F{t,x)) + —{d:^9A{t,F{t,x))Y, 

mi 

dtaB{t,x,y) = dt9B{t,F{t,x),y) + —dx9B{t,F{t,x),y)dx9A{t,F{t,x)) 

mi 


(11.43) 

(11.44) 

(11.45) 

(11.46) 

(11.47) 


and 


dtR{t,x,y) = 


dtTZ(t, F(t, x),y) -\ - dxTZ(t, F(t, x),y)dx9A{t, F{t, x)) 

mi 


2mi 


dtR{t, F{t, x), y)dl9A{t, Fit, x)) 


y^dxF{t,x). 


(11.48) 


Given a solution 

{9a = 9A{t,x),9B = 9B{t,x,y),R = R{t,x,y)) 
of PL35|) - (I1L37|) . we see from (I1L46I) - (I1L48I) that the pair 

[uB = aB{t,x,y),Tl = Tl{t,x,y)) 

satisfies the equations 

dtTZ = —^dyttBdyK - ^OIubR-, 


m2 ^ ^ 

1 x2 1 din 


dtaB = - {dyaB) +„ „ 

m 2 2 m 2 rc 


-V 
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where 


V{t,x,y) := V{F{t,x),y), 

which we recognize as the polar representation of Schrodinger’s equation. This establishes that the wave function ip 
defined by 

ip{t,x,y) := e''°‘^^*'^’^'^TZ{t,x,y) (11.49) 


satisfies the time-dependent Schrbdinger equation 

idtip =--T^d^ip+ V. (11.50) 

2m2 

The above equation shows that the y variable is best interpreted as the configuration variable of a 1-particle 
quantum subsystem in which the motion of the classical particle appears as a time dependent parameter in the 
potential V{t,x,y) = V{F{t,x),y). If V is chosen independent of x, then V = V{y) and (III.501) reduces to precisely 
a 1-particle Schrodinger equation. 

We can further simplify the equation (111.501) by replacing the flow F{t, x) with the Hamiltonian flow associated with 
the Hamilton-Jacobi equation (IH.36p . To see this, consider the Hamiltonian 

Hiix,p) = + U{x) 

and let 


F{t,x,p) = {X{t,x,p),P{t,x,p)) 
denote the correponding Hamiltonian flow, i.e. T satisfies 

r)]-f 

dtX{t,x,p) =-^{F{t,x,p)), (11.51) 

dtP{t,x,p) = — ■^{F{t,x,p)), (11.52) 

F{0,x,p) = ix,p). (11.53) 

Then from the method of characteristics, see [HI, §3.2 & 3.3], it follows that 

F{t, x) = X{t, X, dMO, x)) , (11.54) 

and that the solution ffs to the Hamilton-Jacobi equation (III.361) is given in terms of the Hamiltonian flow by the 
formula 


0Ait,x) - 9A{0,G{t,x)) 


2toi 


\P{T,x,dx6A{0,x))\^ -U{F{t,x)) ) dr 


x—G{t,x) 


(11.55) 


where G{t,x) is for constant t the inverse of F{t,x), i.e. G{t,F{t,x)) = x. 

From the above calculations, it is clear that solutions to the classical-quantum limit system can be generated from 
solutions to the time-dependent 1-particle Schrbdinger equation (III.501) and the classical Hamiltonian flow equations 
(IH.51|) - (|H.53|) via the transformations (III.381) . (III.441) . (111.451) . (III.491) . and (III.541) . and the formula (111.551) . In this way, 
we can view equations (III.501) and (IH.51I) - (IH.53I) as an equivalent formulation of the classical-quantum limit equations 
at least on time intervals for which F{t,x) is invertible, that is before the onset of caustics. Moreover, since (III.501) is 
linear, it is the time of existence for the classical Hamiltonian flow and the formation time for caustics that limit our 
ability to solve the classical-quantum limit equations. 


III. THE ORIGIN OF NONLOCAL SIGNALLING IN THE HALL-REGINATTO HYBRID SCHEME 

To understand the origin of the nonlocal signalling in the Hall-Reginatto hybrid scheme, we consider the case where 
the classical-quantum limit equations (IH.29I) - (IH.30I) arise from the limit of two commuting 1-particle Hamiltonians 
by assuming that the potential H is a function of y only, that is 


V = V{y). 
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We then write the classical-quantum limit equations as 


where 


and 


dt 



Xc{p,e) + XQ{R,e) 


Xc{p,e) = 

( 



\ 2mi 

{da,Aef + 2da,A9d^I3e 



Xq{p,0) = 




B 


m 2 

2 


-ikM + 


1 Dj-JP 

2m2 ^/p 


-V 


(IILl) 


(IIL2) 


(IIL3) 


Since equation (IIII.lIl arises from the limit of a 2-particle quantum system whose evolution is governed by the sum of 
the two commuting 1-particle Hamiltonian operators 


Bi 


1 


2 mi 


-[/ and Ho = — 




and the e \ 0 limit of the 1-particle Hamiltonian operators Hi and H 2 , in our variables, are given by the two vector 
fields (IIH.2I) and (IHI.3I) . respectively, it follows that the flows of these vector fields must commute, and moreover, 
from the discussion above, that Xc and Xq generate the evolution of independent classical and quantum subsystems, 
respectively. Thus, if we let , F^^, and denote the flows of the vector fields Xc, Xq, and Xq + Xq, 

respectively, then F^^ and F^'^ must commute, i.e. 


and is given by 


Ft^^ o Ff« = F^^ o F^'^, 




(HL4) 


(HL5) 


which allows us to write solutions of pil.lll as 

The two properties (IHI.4I) and (|HI.5|) reflect nothing more than the independence of the evolution of the classical and 
quantum subsystems. 

Next, since the quantum system satisfies the no non-local signalling property with respect to position measurements 
of the 1-particle subsystems by virtue of the evolution being generate by the sum Hi + H 2 of commuting Hamiltonian 
operators, it follows automatically that the limit system must also satisfy this property since the position measurements 
are represented by e independent maps in our formalism. In contrast to our equations (IIII.lIl for non-interacting 
classical-quantum systems, those of Hall and Renginatto are given by [l5l| 

dt(^^ =XHRip,0) (H 1 . 6 ) 


where 


Xhr{p,0) = 


These equations are known 


to suffer from a non-local signalling effect with respect to position measurements. 
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Writing the Hall-Reginatto equations as 


where 


dt Q = 0) + Xq{p, 9) + Xi{p, 9) 


Xi{Py := Xhr — Xc — Xq = 


,r^Ad^B9f + A[-^^{dy9) 


1 

2 m 2 VP 


we see that Xj represents a non-local interaction term that alters the evolution generated by the vector field Xc + Xq. 
From this, we conclude that the non-local interaction term Xj is responsible for the non-local signaling effect in the 
Hall-Reginatto equations (IIIL6I) . It is also clear from the above calculations that we can obtain any of Hall-Reginatto 
classical-quantum hybrid equations from the articles [13, [lEIHl by adding a suitable interaction term to Xc +Xq. 


IV. THE FIRST ORDER CORRECTION 


As discussed above, the classical-quantum limit equations derived in Section fll HI may be viewed as the zeroth order 
equations in an e-expansion. For the purpose of improving the approximation and also determining the backreaction 
of the second subsystem on the first, it is necessary to include the effects of higher order corrections. In this section, 
we derive the equations satisfied by the first order correction. 

We begin the derivation by decomposing 9^, and 9% as 


R'^ = R + e/r*^. 


9a = 9a + and 9% = 9b 


(IV.l) 


where and satisfy Av'^ = and Auj’^ = 0, respectively, and, as above, {R,9a,9b) satisfy the classical-quantum 
limit equations (IIL31I) - (IH.33I) . Substituting pV.ll) into the equations (IH.24I) - (IH.26I) . we find, with the help of PL31I) - 
(IH.33|) . that the triple ,u!q) satisfies the following equations: 


5t/i" = -^(dx9AdxP'' + {dxv^ + dxOB + {d^R + {plOAp" 


+ {dlv^ + 919 b + (i? -f ep^)'^ - + edyUJ^)dyp^ - + edluj^)p\ (IV.2) 


dtiy^ = - 


dtuj^ = B 


m 2 " ^ - 2 to 2 

■^—(28x9 aOx^^ + edOxty")"^ + A{dx9B + edxUJ^)^)^ 

\ e (9xR + e^d‘^p\ ^ (pj Q I I ^ PdyR + edyP 

1 f d^R + ed^p'^ 


FA 

1 

2m 


“ i^dx9 aQx^'^ + ‘2,dxi'dx9 B + 2edxio'^dxV + {8x9 b + £8x0j'')^^ + 


R + ep^ 


2m 


— {28y9B8yUj'' + e{8yU}’'f) + 


Y 8yR + e8yp^ 


y R ~\- 

R )\\ 


(IV.3) 


(IV.4) 


We observe that the equations (IH.31F(IIL33ll and (IIV.2ll - (IIV.4l) for the variables {R,9a,9b, p'^,i^'^,uj'^) are equivalent 
to the Schrodinger equation (IlLl) via the transformation 
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Taking the formal e \ 0 limit of the equations (IIV.2I) - (IIV.4|) . we obtain the equations 
dtfJ- = + {dxV + dxOB)dxRj 

— + {dlv + dl9B)R] - —^y0B^y^J. - 

ii \ / m2 Zm2 


2mi 


dtv = - dx0AdxV + A 

mi 


„ 1 dyR 

-{dy0B? +-^ 


dtUj = B 


2m2'“^ ' 2m2 R 

--^^(2dx0AdxiA! + 2dxvdx0B + {dx0B)^^ 


1 d^R 1 ^ ^ ^ 

2toi R 7712^^ ^ 2m2 \ R R^ J 


_ 0d^R \ 


(IV.5) 

(IV.6) 


(IV.7) 


These equations define the first order correction to the limit equations, and we can expect, if 0^^, 9^) is a solution 
of the full quantum mechanical system, then 

= i? + e/x + O(e^), + O(e^) and 0% = 0b + + 0{e^) 

on some fixed time interval [0,T] where {R,9a,0b) solves (IIL31|) - (III.33I) and (^, xx, w) solves (IIV.5l) - (IIV.7ll . Higher 
order corrections can be determined in the usual manner. 

We emphasize that the system consisting of the equations (III.31|) - (III.33I) and (IIV.5I1 - (IIV.7I) represents a fully con¬ 
sistent first order approximation to the Schrddinger equation that should be accurate for modeling the interaction 
between particles of disparate masses in the regime where the large mass particles are effectively classical. We expect 
these equations should be of use in computational chemistry and also the theory of measurement where interactions 
between quantum and classical systems play a distinguished role. 


V. DISCUSSION AND OUTLOOK 

In this article, we have introduced a new approach to analyzing the interaction between classical and quantum 
systems that is based on taking the classical-quantum limit of multi-particle Schrddinger equations. The classical- 
quantum limit equations obtained by this procedure describe the interaction between classical and quantum systems 
and possess many desirable properties that are inherited from the multi-particle quantum systems. 

We considered two applications of the classical-quantum limit equations with the first being to identify the source 
of the non-local signalling that is present in the “non-interacting” classical-quantum hybrid equations of Hall and 
Reginatto. More specifically, we showed that the non-local signalling in the Hall-Reginatto scheme is due to a non¬ 
local correction to our classical-quantum limit equations, which we contend should be taken as the “true” form for 
the equations that govern non-interacting classical-quantum systems in the Hall-Reginatto hybrid scheme. 

In the second application, we viewed the classical-quantum limit equations as the equations satisfied by the zeroth 
order expansion in e of a fully quantum solution. We then, in the standard fashion, derived the system of equations 
satisfied by the first order correction. Together, the limit equations and their first order correction determine a fully 
consistent first order approximation to the multi-particle Schrddinger equation. We expect that these equations will be 
useful for modeling the interaction between particles of disparate masses in the regime where the large mass particles 
are effectively classical. 

The results derived in this article can be taken in many directions. From the point of view of concrete applications, 
the most important is to develop analytical and numerical techniques to solve the limit equations (III.31I) - (IH.3311 along 
with the first order corrections pV.5p - (IIV.7l) in order to determine their suitability for modelling systems that consist 
of nearly classical subsystems interacting with fully quantum ones. Important examples include molecular reactions 
and the interaction of measuring devices with quantum systems. We are currently working on this and hope to have 
progress to report on in the near future. 
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Appendix A: Equivalence of the classical-quantum limits 

As discussed in Remark FlI. 2 [ there is a classical-quantum limit associated to each operator satisfying (III.6I) - (IIL8I) . 
In this appendix, we show that these classical-quantum limits are equivalent by showing that there is an invertible 
transformation between the limits associated to different operators. An immediate consequence of this is that all 
choices of the operator A lead to physically equivalent classical-quantum limits. 


1. Transformation formula 

To establish the transformation formula and its invertibility, we suppose that A and A are two operators that satisfy 
the properties (III.61) - (III.81) . and we let R = 11 — A and S = I — A denote the complementary projection operators. 
Following the prescription (III.141) . we define 

r := ele^ + se^ 
r ;= eie'^ -f 50% 

where 0*^ and are as defined previously by the polar decomposition (III.131) . These two definitions lead to two 
classical-quantum limits defined by 


(0,i?) := lim(6»%i?%, (A.l) 

€\0 

(0,i?) := lim(d%5%. (A.2) 

€\0 

From the relation pi.l5|) . we observe that 

iir -b 56»" = 0*^ = -A0^ + 5r. (A.3) 

e € 

Applying A to this expression yields 

-AAe^ + (A- AA)e^ = -Ae^ (A.4) 

e ^ ^ e 

where in deriving this we used the property A = A. By assumption, A satisfies dyA = 0, and so, for any function 
^ = (j) for some function 0 = 4'{x). But since A((/)) = (j) for all functions (j) = 4>{x) by assumption, it 

follows that 

AA = A, (A.5) 

and from similar considerations, that 

AA = A. (A.6) 

Using (IA.5I) . we obtain from (IA.4I) that 

ir-be(A-A)6i'^ = A0^ (A.7) 

Taking the limit e 0 gives 

Ae = M. (A.8) 

Substituting (EJl) into (Q gives 

which, in turn, yields 

50 = 50-b (A-A)0. 
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by taking the limit e \ 0. Adding this to (IA. 8 I) then yields the relation 


9 = TO 

(A.9) 

where the operator T is defined by 


f = I +A-A. 

(A.IO) 

The operator T is invertible with inverse given by 


f-^ = 1 +A-A. 

(A.ll) 


To see this, we observe that 

f"^f=(]I +A-A){1 +A-A), 

= 1 + A - A + A - A + AA - ^ + AA, 

= I , (A.12) 

where in deriving the last equality we used the identities = A, A = A, AA = A and AA = A. By a similar 
calculation, we also find that 


ff-^ = 1 . (A.13) 

Together, (IA.12I1 and (jA.1311 show that T is invertible with inverse given by (lA.llll . 

Equation (IA.9I1 with T given by (IA.10|) defines the relationship between the two classical-quantum limits (lA.ll) and 
(IA.2L and establishes the existence of an invertible transformation relating the two limits. As discussed above, this 
establishes the physical equivalence of the limits. 


2. Examples of other operators satisfying (III.6ll - (III.8ll 

There are many operators A that satisfy the properties (III. 6 I) - (III. 8 I) beyond the those defined by (III.51) . A more 
general class is defined by 


H4‘)(x,y) = j a{z)(j){x,z)dz, 

where ct{y) is any function satisfying 

J a{y)dy = 1 . 

We observe that we can recover the operators dehned by (HlKI) by setting 


a{y) = 


{b — a) ^ if y G (a, b) 
0 otherwise 


(A.14) 


(A.15) 


(A.16) 


It is also clear that by choosing functions a{y) satisfying (IA.15I) that are distinct from the step-functions (IA.16|) will 
yield operators that are not of the type (III.51) . but nonetheless satisfy properties pi. 6 |) - (III. 8 l) . 

Another class of operators distinct from (I A. 141) is defined by 


A{(j))[x,y) = (j){x,a) 


(A.17) 


where, here, a S M is an arbitrary, but fixed constant. From the Fundamental Theorem of Calculus, which holds for 
regular enough functions (j), we see that 

d 1 / \ 1 

<f>{x,a) = — / (j){x,z)dz = lim 7 -( / (j){x,z)dz— / (j){x, z) dz] = \im - - / (j){x,z)dz. 

dy Ja y=a b- a\J^ J b\,a b - a 






































15 


From this it is clear that the class (IA.17|) of operators is a limiting case of the operators defined by (|II.5p or equivalently 
(IA.14|) with a defined by (|A.16p . 

It is worth noting at this point that although all operators A satisfying the properties (iim-dim are equivalent from 
the point of view of formal limits, the mathematical properties of the operators can be quite different. For example, the 
operators A defined using averaging procedure (I A. 141) are well-defined, assuming appropriate integrability properties 
for the function a, on functions that are much rougher compared to those on which the operators defined via (|A.17I) 
can act. Considerations like this will likely become important in any attempt to go beyond formal considerations 
and rigorously establish the existence of the classical-quantum limit introduced in this article. We expect that only 
particular choices of the operators A will be suitable for rigorous arguments. 
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